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Abstract. - We investigate the role of discrete-particle noise in interface-controlled Ostwald 
ripening. We introduce the noise within the framework of the Becker-Doring equations, and 
employ both Monte Carlo simulations and direct numerical solution of these equations. We 
find that the noise drives the system towards a unique scaling regime describable by a limiting 
solution of a classical continuum theory due to Lifshitz, Slyozov and Wagner. The convergence 
towards the scaling solution is extremely slow, and we report a systematic deviation between 
the observed small correction to scaling and a theoretical prediction of this quantity. 



Introduction. Ostwald ripening (OR) is a generic coarsening process which occurs 
in a late stage of phase separation, when the domains (or clusters) of the minority phase 
compete for monomers. As a result, the larger clusters grow at the expense of the smaller 
ones. Following the pioneering work of Lifshitz and Slyozov [1] and Wagner [2], experimental 
and theoretical investigations of OR have focused on the dynamic scaling properties of the 
probability distribution function (PDF) of cluster sizes, and of its moments. These properties 
are determined by the kinetics of the monomer transport. Simple limits of this kinetics are 
observed when the monomer transport is controlled either by diffusion of the monomers in the 
bulk, or by the processes of attachment and detachment of the monomers at cluster interfaces. 
The latter limit is called interface-controlled. In each of the two limits continuum mean-field 
theories have been formulated: for diffusion-controlled OR by Lifshitz and Slyozov (LS) [1], 
and for interface-controlled OR by Wagner [2]. The two models are often united under the 
name of the LSW model. The LSW model admits a family of self-similar solutions for the 
PDF of cluster sizes. However, the problem of selection of the correct self-similar solution is 
non-trivial. All of the LSW scaling functions have compact support and can be parameterized 
by the value, A, of the logarithmic derivative at the edge of the support; — 1 < A < oo. LS [1] 
argued that the selected PDF for the case of initial conditions with a long tail is the limiting 
PDF, corresponding to A = oo. The scaled PDF should approach this function as t — ► oo. On 
the other hand, if the initial data for the PDF has compact support the selected self-similar 
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PDF is determined by the behavior of the initial data near the edge of its support [3-5] . If 
the initial PDF has logarithmic derivative Ao at the edge of support, the selected self-similar 
solution is the one with the same value of the logarithmic derivative at the edge of its support: 
A = Ao [3,4]. If the logarithmic derivative of the PDF at the edge of support at t = does 
not exist, the PDF does not approach any self-similar solution [5]. 

These results are in apparent contradiction to experimental results which appear to show 
strong selection, i.e. selection insensitive to initial conditions. To find strong selection we 
must consider dynamics beyond that of the classical LSW model. One possibility is to account 
for discreteness of atoms in any real system. This direction was explored by Velazquez [6] 
and Meerson [7], who employed the Becker-Doring (BD) equations [8], properly modified to 
account for conservation of the total number of atoms [9, 10]. In the limit of s S> 1, where s 
is the average number of atoms in a cluster, the BD-equations reduce to the LSW model [9] . 
Taking into account the next order term in 1/s, one arrives at a Fokker-Planck (FP) equation 
for the cluster size PDF [6,7], see below. The drift term of the FP-equation coincides with 
that of the LSW model, while the small diffusion term comes from discrete-particle noise. 
The diffusion term produces a tail in the PDF, even if the initial data has compact support 
at t = 0. According to Refs. [6,7] (see also Ref. [11]), this tail drives the system towards the 
limiting self-similar solution corresponding to A = oo. 

These arguments assume, however, that the FP-equation is a faithful long-time description 
of the BD-equations. Though natural, this assumption is not obviously correct. There are 
many examples when the FP equation misses important aspects of discrete systems [12]. In 
this work we investigate the role of discrete-particle noise by dealing directly with the BD 
equations, without making the FP approximation. We focus on interface-controlled kinet- 
ics. This choice is motivated by experimental findings which showed the importance of this 
(sometimes overlooked) limit in a variety of environments, such as the coarsening of two- 
dimensional islands on Si(001) [13], coarsening of granular clusters in electrostatically driven 
granular powders [14], etc. There is an important additional motivation. Both the LSW 
model and the BD-equations neglect spatial correlations. As a result, the quantitative va- 
lidity of the LSW model is limited, in the case of the diffusion-controlled OR, to extremely 
small area fractions [15]. By contrast, spatial correlations in interface-controlled OR are much 
weaker [16]. Therefore, for the same value of area fraction, the BD-equations are more accurate 
in describing interface-controlled OR, than diffusion-controlled OR [17]. 

BD-equations and the LSW model. - Let s be the number of atoms in a cluster, and 
N s (t) be the number of clusters of size s. The BD-equations [8-10] are master equations for 
the populations of clusters, s > 2: 

N s =N 1 {K S _ 1 N S _ 1 -K S N S )- — + ^±i, (1) 

T~s T s+ l 



and monomers: 



m = -2K x Nl -iV^W + ^+y^. (2) 

s>2 z s>3 

Here K s — K\s p is the rate of attachment of monomers to the cluster of size s, and r s = as q is 
the inverse rate of detachment of monomers from the cluster of size s. One can always choose 
a = 1: this corresponds to a rescaled time t = t/a (the tilde will be omitted in the following), 
and a rescaled attachment rate coefficient a = K x a. Equations and @ preserve the total 
number of atoms N: 

N 1+ J2 sN * = N - (3) 
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To set the stage for our analysis, let us briefly review the predictions of the underlying 
continuum theories. When the dynamics © and © reach the stage of OR, one can proceed to 
the limit of s 3> 1 and treat s as a continuum variable (except for the monomers, s = 1, which 
should be taken care of separately). Then, by a truncated Taylor expansion, one obtains the 
FP-equation [6,7] 

dn s d ,,_ . 1 d 2 ,„ 

!H + 0- S ^ n ^2d^^ n ^ (4) 

where 

V s (t) = An x s p - s- q and D s {t) = An lS p + s~ q (5) 

are the drift velocity and the diffusion coefficient in s-space, n s (t) = N s (t)/N and A = K\aN. 

The LSW model neglects the diffusion term in Eq. an d deals with the continuity 
equation 

^ + £[(W---«H]=o I (6) 

combined with the conservation law 

POO 

ni+ / sn s (t) ds = 1 , (7) 
Jo 

see Ref. [7] for details. For interface-controlled OR one obtains p — 1/2 and q — in two 
dimensions, and p = 2/3 and g = — 1/3 in three dimensions [7]. In the following we will focus 
on the two-dimensional case. At late times the contribution of the monomers to the total 
number of atoms in the system becomes negligible, and Eq. Q becomes 



sn s (t)ds~l. (8) 



It follows from Eqs. JBJ and © that 



•s 



= {&*) = J °- (9) 

Jo n s ds 



where s c = (Ani) -2 is the time-dependent critical size of the clusters; that is clusters with 
s > s c grow, while clusters with s < s c shrink. Equations © and JSJ admit a family of 
self-similar solutions: 

1 , / s \ (3 



,-2 



Ul = MV* ' (10) 



parameterized by (3 — const. According to Ref. [6,7], the discrete-particle noise selects the 
limiting solution with (3 — 2. The corresponding scaling function [2,18] 

if x > 4 , 

has an infinite number of vanishing derivatives at the edge of its support x max = 4 and 
A = oo. In Eq. (JTTI C = e~ 2 + 2 Ei(— 2) = 26.6423 . . ., where Ei is the exponential integral 
function [19]. 
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Fig. 1 - The number density of monomers n\ (a) and the average cluster size s a (b) versus time. The 
circles are averages over 10 3 realizations, obtained in Monte Carlo simulations with N = 3, 001 x 10 6 
atoms. The parameter a = 10 -4 . The initial conditions are: iVi = 10 3 , jV 50 = 6 x 10 4 , the rest of iV s 
is zero. The red lines show the theoretical prediction from Eq. 1101 with f3 = 2 (a) and a linear fit of 
the data (b). 



Monte- Carlo simulations. - Direct Monte-Carlo simulations of interface-controlled OR 
are inefficient, because most of the computation time is spent to resolve the fast monomer 
diffusion in the bulk, while the kinetic bottleneck here is the slow attachment and detachment 
of monomers at the clusters interface [20]. To accelerate the simulations we assume that 
the monomer transport in the bulk is instantaneous. Therefore, we deal with a collection of 
clusters, s > 2, and a pool of monomers, without taking care of the spatial distribution of 
any of them. Essentially, this corresponds to a stochastic simulation of the BD equations. 
We start with a collection of clusters in some initial condition, and some initial number of 
monomers, Ni. We repeatedly choose a cluster at random and let it either grow, by absorbing 
monomers with rate aiVis 1 / 2 , or shrink by emitting monomers with unit rate. 

This is achieved by rejection Monte Carlo: we calculate the rate for the largest cluster, 

1 1% 

s = s maXl to either grow or shrink, f ma x = aNiSmax + 1, and the corresponding quantity for 
the cluster at hand, / = cviVis 1 / 2 + 1. Then we choose a random number u between and 
fmax- If u > f we do nothing. If 1 < u < / we make the cluster grow by increasing s by 1, 
and decreasing Ni by 1. If u < 1 we shrink the cluster, that is decrease s by 1 and increase 
Ni by 1. We need to adjoin a special rule for clusters of size 2. If they shrink we add two 
monomers to Ni and remove the cluster in question from the list. The time was advanced, at 
each step, by St = l/(N c if max ), where N c i is the current total number of clusters (excluding 
monomers) . 

Typical results of these Monte Carlo simulations are shown in Figs. 1 and 2. Figure la 
depicts the number density of monomers n\ versus time, and the theoretical prediction for 
ni from Eq. 1|10(1 with (3 = 2. Good agreement is observed. The average cluster size, which 
we define as s a — (s 1 / 2 ) 2 is plotted versus time in Fig. lb. Linear growth is observed as 
expected. The observed slope 0.268, however, is slightly higher than the theoretical value 
0.25. The scaled PDFs of the cluster sizes, s 2 (i) n s (t), are plotted versus x — s/s c (t) at three 
different times in Fig. 2. Only times up to t — 3, 200 are used, because at later times the 
measured PDFs become too noisy. The three scaled PDFs in Fig. 2 show a good collapse, 
but the scaling function is slightly different from the theoretical prediction Ijlip. Where does 
the deviation come from? The LSW-problem notoriously has very slow convergence to the 
limiting self-similar solution [1,21]. Therefore, one can attempt to attribute the observed 
deviation to very slow convergence, masked by fluctuations. 
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Fig. 2 - The scaled PDFs of the cluster sizes at times 400 (black), 800 (blue) and 3200 (green) as 
obtained in Monte Carlo simulations. Because of the relatively large noise at late times the black and 
blue lines are masked by the green line. The parameters are the same is in Fig. 1. The red dashed 
line is the theoretical scaling function J I I i . 

Numerical solution of the BD- equations. To test this interpretation and reach later 
times, we solved the BD-equations numerically. Equations Q and (0) for 1 < s < 10 5 were 
solved using a fourth order Runge-Kutta algorithm. The conservation law was used for 
accuracy control. The time step chosen was At = min [0.007 /(aN\), 0.125]. This choice 
resulted in excellent mass conservation and enabled us to probe very long times. For the 
results reported below only 10~ 7 of a single particle was lost by t = 10 5 . 

The results are shown in Figs. 3 and 4. The number density of monomers n\ versus 
time (Fig. 3a) agrees very well with the theoretical prediction. The average cluster size 
s a = (s 1 / 2 ) 2 shows linear growth with time, see Fig. 3b. The slope 0.262, obtained on the 
interval 10 3 < t < 10 5 , is still slightly higher than the theoretical value 0.25. The scaled PDF 
(Fig. 4a), though still different from the theoretical scaling function apparently slowly 
approaches it. 

Logarithmic corrections to scaling. LS [1] investigated, in the context of diffusion- 
controlled OR in three dimensions, the convergence towards the limiting self-similar solution. 




Fig. 3 - The number density of monomers rn (a) and the average cluster size s a (b) versus time, 
obtained by solving Eqs. Q and @ numerically for N — 3 x 10 6 . The parameter a — 10 -4 . The 
initial conditions are: Ni = 1.5 x 10 3 , N 50 = 59 9 70, the rest of N s is zero. The thin red lines show a 
linear fit of the data (a) and the theoretical prediction from Eq. (I1UB with /3 = 2 (b). 
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Fig. 4 - The solid lines in (a) show the scaled PDFs at times 1350 (black), 10350 (blue) and 85350 
(green), obtained by solving Eqs. Q and numerically. The parameters are the same is in Fig. 3. 
The red dashed line is the theoretical scaling function 1111 . The black dashed line in (b) shows A(t) 
(see text), while the red solid line shows ln~ 2 t, see Eq. p5^ . 



Here we will employ their argument (see also Ref. [6]) in the problem of interface-controlled 
OR in two dimensions, and compare it with our simulations. 
The characteristics of Eq. JBJ are described by the equation 

a = (s/sc) 1/2 - 1 . (12) 

In the rescaled variable x — s/s c Eq. (|12f> becomes 

s c ± = -s c x + x 1/2 - 1 . (13) 

The limiting self-similar solution requires that s c (t) — > i/4 as t — > oo, so that the right hand 
side of Eq. l(T3j l becomes a perfect square: tx = —(^fx — 2) 2 . As shown by LS, this is 
the only possibility to have a non-diverging normalization integral (JBJ), when the full PDF 
has a tail. Following LS, we are looking for the leading correction in the following form: 
s c (t) = it/4) [1 + e 2 (t)], where e(t) <C 1. Consider a small region of x around the "blocking 
point" x = 4 (the edge of support of the similarity solution). Let y(t) = [x(t) — 4] /s(t). 
Equations l|13|) becomes, in the leading order, 

±dy = _ 1 y±(l)_l_ (14) 

4edr 4dT\e) 64' 1 ' 

where r = Int. By the same normalization argument, the right hand side should become a 
perfect square as t — > oo, which yields e = 1/r = 1/lni. Therefore, we obtain 

.„ = <,VV = -„(*) = | + + ■ (15) 

The logarithmic correction ln~ 2 1 implies that the apparent dynamic exponent would be 
slightly larger than 0.25 and decreasing in time, as indeed observed in our simulations. How- 
ever, a more detailed comparison of Eq. H15J) with our simulation data is dissapointing. Figure 
4b shows the quantity 8 — 4s c (t)/t—l versus t. Though 6 does go down with time (apparently, 
logarithmically slowly), it clearly disagrees, at these times, with the prediction of Eq. |TH]l. 
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Summary. - We have investigated interface-controlled OR in the framework of conserved 
BD-cquations. We performed Monte-Carlo simulations of the system, and also solved the BD 
equations numerically. We observed that discrete-particle noise drives the system towards 
a limiting self-similarity solution described by the LSW theory. However, the convergence 
towards the scaling regime is extremely slow. Furthermore, there is a clear disagreement 
between the observed small correction to scaling and a theoretical prediction of this quantity, 
obtained in the spirit of the LS theory. At present we cannot pinpoint the reason for the 
disagreement in the subleading order of the theory. Clarifying this important issue should be 
the next step of theory. 
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